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Abstract. Symmetries and conservation laws are studied for two classes of physically and 
analytically interesting radial wave equations with power nonlinearities in multi-dimensions. 
The results consist of two main classifications: all symmetries of point type and all conserva- 
tion laws of a general energy-momentum type are explicitly determined, including those such 
£jrj as dilations, inversions, similarity energies and conformal energies that exist only for special 

^ 1 powers or dimensions. In particular, all variational cases (when a Lagrangian formulation 

exists) and non- variational cases (when no Lagrangian exists) for these wave equations are 
considered. As main results, the classification yields generalized energies and radial momenta 
in certain non-variational cases, which are shown to arise from a new type of Morawetz di- 
lation identity that produces conservation laws for each of the two classes of wave equations 
in a different way than Noether's theorem. 

1. Introduction 

Symmetries and conservation laws are important tools in the study of global analysis of 
nonlinear wave equations 

u a = g(x, u, Vm) ■ (y 2 u) + f(x, u, Vu) (1) 

for u(t, x) in n > 1 spatial dimensions (where x denotes Cartesian coordinates). 

As pointed out in [161 E| > conservation laws such as energy provide basic conserved quan- 
tities used in obtaining estimates on \u\ or |Vw| for classical solutions, and also in defining 
suitable norms for weak solutions. Of considerable interest are extra conservation laws such 
as conformal energies that can appear in the case of power nonlinearities \u\ p or derivative 
nonlinearities \Wu\ p in g and / for special powers p depending on the dimension n. Symme- 
tries, in contrast, lead to exact group- invariant solutions and play a role in defining invariant 
Sobolev norms. Scaling symmetries are of special relevance, as the critical nonlinearity power 
for a blow-up is typically singled out by scaling-invariance of a positive energy norm. More- 
over, scaling transformation arguments give a means of relating the behavior of solutions in 
different regimes, for instance, solutions at short times with large initial data can be scaled 
to long times with small initial data when the nonlinearity power is subcritical. 

In previous work by Anco&Ivanova (5], conservation laws and symmetries have been clas- 
sified for semilinear wave equations 

u tt = Au + au v (2) 

with a power nonlinearity in n > 1 spatial dimensions. In the present paper and a sequel, we 
extend those classifications to a variety of physically and analytically interesting quasilinear 
wave equations 

^ = V- ((c + V)Vu) + (b--)u p Au (3a) 
l 



and 

uu = V • ((c+ ^|Vu| p ) Vu) + (6 - -) |Vu| p Au (3b) 

contained in the general class (pD), with free constants a, b, c, where p 7^ is the nonlinearity 
power and n > 1 is the spatial dimension. These wave equations can be written respectively 
in the equivalent forms 

u tt = (c + bu p )/\u + au p - l \Vu\ 2 (4a) 

and 

u tt = (c + b\Vu\ p )Au + |a| Vw| p - 2 (Vw • V) 2 w (4b) 

which are valid for any power p. There has been considerable mathematical interest in the 
analysis (e.g. global existence, uniqueness, regularity, blow up) of solutions to the Cauchy 
problem for translationally- invariant quasilinear wave equations in multi-dimensions [TJ [T5] , 
particularly for radially symmetric initial data. The two classes of wave equations fj4aj) - fl4bl) 
studied here provide the simplest such examples having homogeneous power nonlinearities. 

To begin, we will consider the radial reduction of the wave equations f|4"al) and (14bl) for 
u(t,r), where r = \x\ is the radial coordinate in R n . This reduction yields 

Th — 1 

Utt = (c + bu p ) (u rr H u r j + au p ~ l u 2 r (5a) 

and 

Utt = (c + (b + a)u p ) (u rr H u r ) — a u p+1 (5b) 

respectively, which are radial wave equations belonging to the general class 

u tt = g(r,u,u r )u rr + f(r,u,u r ) (6) 



where 



71 — 

g = c + bu p , f = (c + bu p )u r + au p - l u 2 r (7a) 



for equation fl5a|) . and where 



g = c + (b + a)u p , f = (cu r + bu P+1 ) (7b) 

r 

for equation fl5b|) . We note that the class (jTajl is quasilinear hyperbolic iff bp ^ or is 
semilinear hyperbolic iff bp = 0, a ^ and c + b ^ 0; and the class (I7bj) is quasilinear 
hyperbolic iff (b + a)p ^ or is semilinear hyperbolic iff b + a = 0, be ^ and p(p + 1) ^ 0. 

These classes f!7a|) and f jTbj) of radial wave equations describe models for nonlinear ra- 
dial wave propagation in which the wave speed can depend on the wave amplitude or its 
radial gradient. Such models encompass a wide variety of physically interesting examples 
in n = 2 and n = 3 dimensions: for class (!7a|) . radial compressible polytropic fluid flow 
in Lagrangian coordiantes [13J and radial propagation of volt age /current in electromagnetic 
transmission guides [11] : and for class ( !7a|) . radial nonlinear vibrations of membranes and 
radial deformations of hyperelastic materials [10J H2]- The extension of these examples to 
arbitrary n > 1 dimensions is mathematically natural in light of the analytical interest in 
n-dimensional quasilinear wave equations. Consequently, our results in this paper should be 
of direct relevance to investigation of physical models of nonlinear radial wave propagation 
as well as to on-going analytical work on the Cauchy problem for radial quasilinear wave 
equations. 
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There are two important cases in which the radial wave equations (I7a)l and (17b j) have a 
special analytical structure. First, a radial wave equation fl6]) has a n-dimensional divergence 
structure if its derivative terms can be expressed as a total divergence of the form 

(r n - l u t ) t = (r n - l F) r (8) 

for some function F(u,u r ). Such equations ([H]) give rise to kinematic conservation laws 
involving the wave amplitude u and its time derivative Ut, which will be discussed later. 
Second, a radial wave equation ([6]) has a n-dimensional variational structure if it arises as 
the stationary points of a functional 

i*t\ roc r r* 

C= / Lr n ~ l drdt, —=r n -\u tt -gu rr -f) = (9) 

Jt Jo OU 

for some Lagrangian L(u,u t ,u r ). Wave equations with this structure © possess an energy 
conservation law related to time-translation symmetry through Noether's theorem [T4"HT]. as 
will be reviewed in more detail later. 

Wave equation ( |5al) is variational when (and only when) its coefficients a, b, c satisfy the 
relation 2a = pb, where the Lagrangian is given by 

L = \{-u 2 + (c + bu p )u 2 r ). (10a) 

Wave equation (I5b|) is variational when (and only when) its coefficients a, b, c satisfy the 
relation a = pb. The Lagrangian is given by 

L = \(-u 2 t + cu 2 r ) + bh(u r ) (10b) 

where 



h= I u p r +l du r = { p + 2 r ' ^ . (11) 
lnw r , p = —2 

Both wave equations (l5a| and (15b j) are total divergences when (and only when) a = bp, 
with the respective fluxes given by 

F = (c + bu p )u r = {cu + bH(u)) r (12a) 

and 

F = (c + bu p r )u r (12b) 

where 

H = [ u p du = { p+l" 1 r ' ' . (13) 



We note that the conditions for existence of a divergence structure (IB]) and a variational 
structure coincide for wave equation (l5bl) . 

Interestingly, wave equation ( I5b|) has an alternative divergence structure and variational 
structure when c = without any conditions on a and b. In this case 

Tl — 1 

utt = (a + b)u p u rr + b < +1 (14) 

r 
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is a total divergence of the form 




(b + a) \nu r + b(n — 1) lnr, p — —1 



(15) 




I 



oo 



Lr m drdt, L — —\ 



\u\ + / F(r, u r )du r 




(16) 



where 



m = 



b{n - l)(p + l) 
a + 6 



(17) 



Compatibility between variational structures, divergence structures, quasilinearity or semi- 
linearity is summarized by the following result. 

Proposition 1. (i) A radial wave equation (|5a|) has a hyperbolic, quasilinear, variational 
form ([9]) iff '2a = pb =/ ; or a hyperbolic, quasilinear, divergence form ([8]) iff a = pb ^ 0. 
(In particular, these two structures are mutually incompatible) . 

(ii) A radial wave equation ( 15bj) has a hyperbolic, quasilinear, variational/ divergence form 
OH])-© iff a = pb, p(p+l)b ^ 0, or a hyperbolic, quasilinear, alternative variational/ divergence 
form 0151) - flT6"j) iff c = 0, p(a + b) ^ 0. (Moreover, these structures coincide when c = 0, 
a = pb, p(p + 1)6 7^ ; so thus m = n — 1.) 

(hi) yl radial wave equation (I5a|) and (15b j) cannot have a hyperbolic, semilinear, variational 
or divergence form. 

In this paper, for the two classes of radial wave equations (!7al) and (17bl) . the main goals 
of our work will be to determine: 

(i) what symmetries are admitted other than time-translation; 

(ii) what conservation laws are admitted other than energy in the variational case (jH]), and 
kinematic quantities in the divergence case ([H])- 

In section [21 we first recall the definitions of symmetries and conservation laws for the 
general class of radial wave equations §6§. For the case of a radial divergence structure 
([8]) or fTToT) . we discuss the notion of kinematic conservation laws. For the case of a radial 
variational structure (Q or ffl6|) . we give a statement of Noether's theorem and summarize 
the energy conservation law arising from time-translation symmetry. 

In sections [3] and HI we present the classifications of symmetries and conservation laws 
for each class of radial wave equations f !5a|) and 05bl) . For special nonlinearity powers p and 
for special relations among the coefficients a, b, c, both of these wave equations are found 
to possess a variety of interesting symmetries, including scalings, non-rigid radial dilations, 
and a temporal inversion, as shown in section |3j The commutator algebra for all of the 
symmetries is worked out in detail. Variational symmetries are also identified, which yield 
conservation laws through the variational structures for each of the wave equations. As 
shown in section HJ these Noether conservation laws include similarity energies, conformal 
energies, and radial momenta. Interestingly, additional conservation laws are found to exist 
for certain nonlinear powers p and for certain relations among the coefficients a, 6, c. These 
conservation laws describe generalized-energies that do not arise from Noether's theorem. 
We account for them by deriving a variational type of Morawetz radial dilation identity that 
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produces conservation laws from symmetries in a different manner than the usual Noether 
correspondence. This main result is explained in subsections 14. 1 . 2 | and 14.2 .21 for the respective 
wave equations (15a]) and (l5bl) . 

Our classification results will be extended to the original two classes of translationally- 
invariant wave equations ( jjaj) and (lib]) in a subsequent paper. 

2. Preliminaries 

To proceed we state the definitions of symmetries and conservation laws from an ana- 
lytical perspective (see also [7j). A point symmetry of a wave equation (ED is a group of 
transformations given by an infinitesimal generator 

5t = r(t,r,u), Sx = £(t,r,u), 5u = r](t,r,u) (18) 

on the variables t, r, u, under which the wave equation ([6]) is preserved. On solutions u(t, r) 
of equation ([6]), such an infinitesimal transformation ffT8]) is equivalent to a generator 

X = (rt - rut - £u r )d/du (19) 

called the characteristic form of the point symmetry. The scalar functions rj, r, £ are 
determined by a linear equation arising from the Frechet derivative of the wave equation ([6]) 
applied to the function P = rj — ru t — C,u r , 

o = dip - (| + ujj-y - ( ^ + Urr ^ DrP _ gDlP (20) 

holding for all formal solutions u(t,r), namely, with uu (and its r-derivatives) replaced by 
gu rr +f (and its total r-derivatives). (More precisely, symmetries are computed in a jet space 
whose coordinates are defined by t, r, u and all derivatives of u modulo the equation (jSJ) and 
its differential consequences. A point symmetry is then the prolongation of the operator ffl~9]) 
that annihilates the wave equation (jSJ).) 

The definition of a symmetry can be extended to involve first-order derivatives of u by 
considering an infinitesimal generator of the form 

X = P(t,r,u,u t ,u r )d/du (21) 

where P is a scalar function depending on u t , u r , in addition to t, r, u. Such a generator 
corresponds to the infinitesimal transformations [7j 

dP c- 9P 
St = -j—, Sr = -—, (22a) 
ou t OU T 

dP dP 

5u = P-u t —-u r — i (22b) 
ou t ou r 

dP dP 

5u t = — + u t — , (22c) 
at ou 

dP dP 

Su r = — + u r — (22d) 
or ou 

on the variables r, t, u, u t , u r (viewed as coordinates in jet space). A group of transformations 
( 122]) is a contact symmetry of a wave equation <^ if these transformations preserve the 
equation (EJ). 

The set of all infinitesimal point or contact symmetries admitted by a given wave equa- 
tion ([()]) has the structure of a Lie algebra under commutation of the operators X (prolonged 
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to the jet space). For a given (sub)algebra of point or contact symmetries, the corresponding 
group of transformations has a natural action [HJ EE [7J on the set of all solutions u(t,r) of 
the wave equation. A solution u(t, r) is invariant under a one-dimensional (sub)group with a 
generator (|2~T|) if it satisfies the equation P(t, r, u, u t , u r ) = where P = r\ — ru t — £w r in the 
case of a point symmetry. Such solutions of a wave equation (ED are called group invariant. 

A conservation law of a wave equation flS]) is a space-time divergence given by a linear 
combination of the equation and its differential consequences, so that 

D t T + D r X = (23) 

holds for all formal solutions u(t,r), where T and X are scalar functions depending on t, 
r, u, Ut, and r-derivatives of u, u t . The radial integral of the conserved density T formally 
satisfies 

(/ r Tdr = -Xr (24) 



dt 

and so if the radial flux X vanishes both at spatial infinity and at the origin, then 

POO 

C = Tdr = const. (25) 
Jo 

formally yields a conserved quantity for the wave equation (J6]). Conversely, any such con- 
served quantity arises from a conservation law (|23|) . Two conservation laws are equivalent 
if their conserved densities T differ by a total radial derivative D r Q on all formal solutions 
u(t,r), thereby giving the same conserved quantity C up to boundary terms. Correspond- 
ingly, the fluxes X of two equivalent conservation laws differ by a total time derivative —D t Q 
on all formal solutions u(t,r). The set of all conservation laws (up to equivalence) admitted 
by a given wave equation (J6J) forms a vector space, on which there is a natural action [H] by 
the Lie group of all point or contact symmetries of the equation 

Each conservation law (123]) has an equivalent characteristic form given by 

D t T + D r (X — T) = (u t t — gu rr — f)Q (26) 

with 

ST 

■D' T (u tt - gu„ - f) (27) 



E 



i>o 1 1 

where Q is a scalar function that depends on t, r, u, u t , and r-derivatives of u, u t . Such 
functions Q whose product with equation yields a total space-time divergence (and 
hence a conservation law (|23|) on all solutions u(t,r)) are called multipliers. Through the 
characteristic equation (1261) . multipliers and conserved densities (up to equivalence) have a 
one-to-one relation such that Q is the variational derivative of T with respect to u t , 

Q = T~ (28) 
ou t 

which can be inverted to express both T and X (modulo total derivative terms D r Q and 
—DtQ respectively) in terms of Q by means of an explicit integral formula [31 HI [20] based 
on a homotopy integration of the characteristic equation (|26l) . 

All multipliers Q for a wave equation ([6]) are determined by a linear system [HI HI [HI [20] 
that holds for all formal solutions u(t,r). (More precisely, the computation of multipliers 
uses the same jet space as for the computation of symmetries.) The system consists of [7] 
a linear equation given by the adjoint of the Frechet derivative of equation ([6]) applied to 

6 



Q, augmented by additional linear equations formed from the Frechet derivative of Q itself. 
Thus the determination of conservation laws via multipliers is a kind of adjoint problem [2] 
of the determination of symmetries. 

From the variational relation (1281) . conserved densities T up to 1st order correspond to 
multipliers Q of at most the same order, as given by 

Q(i , r! „,^„ r) = ^_WV). (29 ) 

OUt 

Such multipliers are determined by the linear system 

= D tQ - (g + + D r ((§L + u^Q) - Dl( 9 Q) (30) 

= 2* + A*** (31) 
ou ou t ou r 

holding for all solutions u(t,r), namely, with u u (and its r-derivatives) replaced by gu rr + f 
(and its total r-derivatives). For a given multiplier Q, the corresponding conserved density 
T and flux X can be obtained by integration of the relations 

du~ = Q ' du~ = ~ gQ i (32) 

(Jilt ou r 

dT dX n dT dX dT dX 
du r ou t ot or ou ou 

arising from the characteristic equation (|26|) . In particular, an explicit homotopy integration 
formula is given by 

T(t,r,u,u t ,u r ) =Ut / ( Q(t, r, Am, Xu t , Xu r ) — Xu——(t, r, Am, Xu t , Xu r ) )dX 



o 



du 

f 1 ( oW 

+ u r I ( W(t, r, Am, Xu t , Xu r ) + Xu— — (t, r, Xu, Xu t , Xu r ) ) dX 

Jo ^ ou 
(j(r, Xu, Xu r ) — (t, r, Xu, Xu t , Xu r ) 
dQ. A A . . dW 



(34) 



OL) . , . . OW . , _ . \ „ 

+ -—— (t, r, Am, Aw t , AmJ — (t, r, Xu, Xu t , Xu r ) )dX 

Ot or J 

and 

f 1 ( dW \ 

X(t, r, u, u t , u r ) = — u t J [W{t, r, Xu, Xu t , Xu r ) + ^ u ~q — (A r ? ^ u t, Xu r )jdX 

n 1 r\ / \ 

— u r j ( (Qg)(t, r, Xu, Xu t , Xu r ) — Xu — - — (t, r, Xu, Xu t , Xu r ) )dX 
J V dM / 



(35) 

dQ 

- a I (^-^-(t,r, Xu, Xu t , Xu r )g(t,r, Xu, Xu t , Xu r ) 

r > Xu,Xu t ,Xu r ) - ^-(t,r, Xu,Xu t ,Xu r j^dX 



where 

W(t, r, u, u r , u t ) = u t — — (t, r, u, Xu t , Xu r )dX + u r g(r, u, Xu r ) — — (t, r, u, Xu t , Xu r )dX. 
Jo ou r Jo du t 

(36) 

These expressions (|34|) - fl36|) are an improvement over the standard homotopy integral for- 
mula in which both T and X would contain r-derivatives of u t , u r that have been canceled 
here by the addition of suitable total derivative terms. (The proof that T and X satisfy 
relations (I3~2l)- (l33"l) assumes that the scaled characteristic expression 

(Xu tt — g(r, Xu, Xu r )Xu rr — f(r, Xu, Xu r ))Q(t, r, Xu, Xu tl Xu T ) 

as well as the scaled differential form 

dQ dQ 

— — (t, r, u, Xu t} Xu r )Xdu t + g(r, u, Xu r )— — (t, r, u, Xu t , Xu r )Xdu r 

are differentiate with respect to A and vanish at A = 0.) 

If a wave equation §6§ possesses a variational formulation then its multipliers Q 
correspond to variational symmetries defined by the characteristic form 

X = r l - n Qd/du (37) 

where the determining equations on Q can be shown to reduce to conditions equivalent to 
those given by Noether's theorem [I^t [71 [8] for the Lagrangian functional £ to be formally 
invariant (up to boundary terms at spatial infinity and the origin) under the prolongation of 
the generator fl37j) . More precisely, symmetry invariance requires that the Lagrangian itself 
satisfies 

r n " 1 prXL = D t A + D r B (38) 

for some scalar functions A and B depending on t, r, u, u t , and r-derivatives of u, u t , where 
X is a symmetry generator (12 ip in characteristic form (and pr denotes its prolongation). 
Then the variational identity 

r"->XL = + Dt^-^P) +D r (r n - 1 ^p) (39) 

yields the equation 



Su V du+ J \ du 



given by 



a r 

—P = (u tt - gu rr - fY^P = D t T + D r X (40) 
ou 

T = A - Q-^—, X = B-Q?-, Q = r n - 1 P (41) 
ou t ou r 

Thus, each symmetry generator ([2T]) for which the Lagrangian functional is invariant (|38|) 
yields a conservation law in characteristic form (|26|) with the conserved density T and flux X 
given in terms of A, B, P by the Noether relation ( 14~T|) . Conversely, through the variational 
identity (139]) combined with equation (T40|) . each conservation law in characteristic form (126]) 
yields a symmetry generator (l37j) under which the Lagrangian functional is invariant (138|) . 
with A and B determined in terms of Q,T,X from the Noether relation (T4T|) . Moreover, 
since invariance of the Lagrangian functional £ implies that its stationary points d~£/5u = 
are preserved, every variational symmetry (1371) of a Lagrangian for a wave equation ([6]) is 
thus a symmetry of the wave equation itself. 



Note that when the multiplier Q given by a Noether conservation law (j4Tj) is of 1st order 
(|29|) then the corresponding variational symmetry (1371) will be of point type if 

9 2 Q = d 2 Q = _^Q_ = Q (42) 
dutdut du r du r dutdu r 

and otherwise it will be of contact type. 

By inspection, any Lagrangian functional of the form fl9]) is formally invariant under in- 
finitesimal time-translation 

St = 1, Sr = 0, Su = 0. (43) 

The corresponding characteristic generator thereby yields a variational symmetry 

X = r^Qd/du, Q = -r n - x u u (44) 

admitted in the case of a variational wave equation (JSJ). The resulting conservation law (JHD 
describes a conserved energy 

C = - Tdr = {u\ + L(u,u t ,u r ))r n - 1 dr (45) 
Jo Jo 

for all solutions u(t, r) that have zero flux at spatial infinity and at the origin. 

In the case of a wave equation f[T4l) having the alternative variational structure (If 6p . the 
previous statements hold with r m in place of r n ~ l . 

For any wave equation (jSj), whether or not it is variational, a conserved density T and a 
flux X that arise from a multiplier with at most linear dependence on u t 

Q = a(t, r, u, u r )ut + (5 (t, r, u, u r ) (46) 

will be called a generalized- energy type conservation law. For such conservation laws, the 
relations ( I32p - (|3"3"|) show that the form of T and X is given by 

T =~u 2 a(t, r, u, u r ) + u t /3(t, r, u, u r ) + A(t, r, m, w r ) 



i 3 da(t,r,u,u r ) l 2 d/3(t,r,u,u r ) dA(t,r,u,u r ) (47) 

X = ~ ^ ^ ^ ^ Ut d£ + B(t > Mr) 

for some scalar functions A and B. 

If a wave equation (jSJ) possesses a divergence structure OH]), then it can be integrated to 
get 

j2 /-oo 

JLjf ur"- 1 cir = F| 00 (48) 
which formally yields a kinematic quantity given by 

poo 

C(t) = / ur n ~ l dr = Cit + C 2 (49) 



satisfying 

provided that the flux F vanishes at spatial infinity and at the origin. The resulting quantities 

POO 

C\ — I u t r n ~ l dr = const. (51) 
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o 



C 2 = (u- tutY^ 1 dr = const. (52) 
Jo 

are conserved, corresponding to the conservation laws 

T\ = r n - l u t , X x = -r n ^F (53) 

T 2 = r n -\u - tut) , X 2 = r n ~HF (54) 
whose multipliers are given by 

Qi = r n ~\ Q 2 = ~r n -H. (55) 

In the case of a wave equation (IT4T) having the alternative divergence structure (115]) . the 
same kinematic conservation laws hold with r m in place of r"" 1 . 

More generally, for any wave equation (E]) whether or not it has a divergence structure, 
a conserved density T and a flux X arising from a multiplier Q that depends only on the 
coordinates t, r will be called a kinematic type conservation law. Such conservation laws are 
characterized by the form 

9a(t,r) , . 

T = u t a<(t, r) - u — — — (56) 

for the conserved density, where 

Q = a{t,r) (57) 

is the corresponding multiplier. 



3. Symmetry Classification 

For each class of radial wave equations ( I5al) and ( I5bj) we will now find all point symmetries 
( 1371) . In particular we will explicitly determine any point symmetries that exist only for 
special nonlinearity powers p and dimensions n / 1, as well as for special relations among 
the constant coefficients a, b, c in these wave equations, excluding cases where the wave 
equation is linear or non-hyperbolic. (Note we will allow n to have non-integer values. An 
interpretation of the wave equations ([5a]) and (15bj) in such cases is given in section [5j) 

Our results are obtained by solving the determining equation ( 120]) for the functions 
r(t, r, u), £(t, r, u), r](t, r, u) given by the characteristic form of the symmetry generator (IIU]) . 
To keep the list of solutions succinct, we have merged symmetries with similar forms into a 
combined form wherever possible. Remarks on the computation will be provided at the end 
in section 13.31 

We will also work out the commutator structure of these symmetry generators. For this 
purpose it is convenient to use the associated canonical form of the generator defined by 

X = rd/dt + id/dr + r]d/du (58) 

which directly corresponds to an infinitesimal transformation ( 1T8]) on (t,r,u), namely 5t = 

X.t = r(t, r, u) , 5r = Xr = £(t, r, u), 5u = Xu = r](t, r, u). 
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|p6r 


6w + pc 


cp(p- 1) =0, 6^0 


4 





(a + 6 + (1 -p)c)r 3 ~ n 


(2 - n)r 2 ~ n ((6 + (1 - p)c)u + pc) 


cp(p — 1) = 0, 

(a + b + (1 - p)c)(n - 3) 

- op{n/Z - i), 
n^2 


5 





|p6r In r 


(1 + lnr)(fow +pc) 


cp(p - 1) = 0, 

a = (p — l)c — b, b 7^ 0, 

n = 2 


6 


t 2 







c = 0, 
p = -4 


7 








1 


6 = 0, 
p = 1 


8 








t 


6 = 0, 
p = 1 



Table 1. Infinitesimal point symmetries for 1, (6p) 2 + a 2 (c + 6) 2 7^ 0. 



3.1. Point symmetries of u u = (c + bu p )(u rr + (n — l)u r /r) + au p ~ l u 2 . 

Tabled] lists the infinitesimal point symmetries of the radial wave equation (I5al) in all cases 
such that this equation is nonlinear (i.e. a 7^ or pb 7^ 0), hyperbolic (i.e. p6 7^0orc + 67^0 
when pb = 0), and multi-dimensional (i.e. n / 1). These three restrictions are equivalent 
to the inequality (n — l)((bp) 2 + a 2 (c + 6) 2 ) 7^ holding on the nonlinearity power p, the 
dimension n, and the coefficients a, 6, c. 

Symmetries Xi and X 2 respectively generate a time-translation and a space-time dilation 
on (t,r). These are the only symmetries admitted for all allowed values of a, 6, c, p, n. 
Symmetries X 7 and X 8 generate a shift on u, while symmetry X 6 generates a temporal 
inversion on (t,u). Symmetry X 3 generates a scaling on (r,u) when c = 0, or a scaling on 
u when p = 0, and otherwise generates a shift on u combined with a scaling on (r,u) when 
c ^ 0, p ^ 0. Symmetry X4 is a non-rigid generalization of the shift-scaling X3 on (r, u). 
Symmetry X5 is a logarithmic counterpart of X4. 

From Table [1] we see that the point symmetry structure for the wave equation (l5a| is 
richest in four main cases: n = 2, i.e. when the spatial domain is planar; p = 1, i.e. when 
the nonlinearities are quadratic; p = 0, i.e. when the Laplacian term is linear; and c = 0, i.e. 
when the linear Laplacian term is absent. We also see that the only distinguished relation 
among a, 6, c is given by (a+6+(l— p)c){n— 3) = bp(n/2 — l), which reduces to a = (p-l)c-b 
when n = 2, where p and c are restricted by the condition p(p — l)c = 0. When the radial 
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wave equation (!5a| is quasilinear (i.e. p 7^ 0, 6 7^ 0), this relation becomes 

a = ~{p l ~ n ^ +l)b, n^3 (59) 
V n — 3 J 

with p = 1 or c = 0. In contrast, when the radial wave equation (15 aj) is semilinear (i.e. 
p6 = 0, a^O), the distinguished relation among a, 6, c simplifies to 

a = —(6 + c) or n = 3 (60) 

We note that, from Proposition [H the simplified relation ( 159]) in the quasilinear case coincides 
with the condition 2a = pb for existence of a variational structure (Q if p = 2{n — 3) 7^ or 
the condition a = pb for existence of a divergence structure §8§ if p = 2{n — 3)/(4 — n) 7^ 0, 
while the non-existence of both a variational structure and a divergence structure in 
the semilinear case precludes any overlap with the relation (160]) . 

We will next summarize the structure of the algebra generated by all of these point sym- 
metries. 

Theorem 1. For a multi- dimensional, hyperbolic, quasilinear radial wave equation (15a]) . its 
point symmetries in the general case n ^ 1, p ^ 0, 6 7^ comprise a 2-dimensional algebra 
generated by the time-translation X4 and the space-time dilation X 2; with the commutator 
structure 

[X 1 ,X 2 ]=X 1 (61) 

Its additional point symmetries generate larger algebras in the following cases: 
(i) For c = 0, the point symmetries additionally comprise the scaling X 3 ; the non-rigid 
dilation X 4 when a(n — 3) + b(n — 3 — p{n — 2)/2) = 0, n ^ 2; the logarithmic dilation X 5 
when a + b = 0, n = 2; and the temporal inversion X 6 when p = —4. In the subcases n = 3, 
P 0; P ~ 4; or n / 3, a ^ — b(j>{\ — n/2)/{n — 3) + l), p 7^ —A, the scaling X 3 enlarges 
the commutator structure (IBT]) by 

[X 1 ,X 3 ] = [X 2 ,X 3 ] = (62) 

generating a 3-dimensional algebra. In the subcases n ± 2, n ^ 3, a = -b(p(l - n/2)/{n - 
3) + l) , p 7^ — 4, or n = 3, p = 0, the non-rigid dilation X4 enlarges the commutator structure 
flMD-il by 

[X l5 X 4 ] = 0, [X 2 , X 4 ] = (2 - n)X 4 , [X 3 , X 4 ] = p(l - n/2)6X 4 . (63) 

In the subcase n = 2, a = —b, p 7^ —A, the logarithmic dilation X 5 enlarges the commutator 
structure (IBT]) - fl6"2"]) fry 

[X 1; X 5 ] = 0, [X 2 , Xb] = X 3 , [X 3 , X 5 ] = pb/2%. (64) 

In the subcase n^2, n^3, — b(j>{\ — n/2) /{n — 3) + l), p — — 4 ; the temporal inversion 
X 6 enlarges the commutator structure (IBT]) - (162]) fry 

[X 1; X 6 ] = 2X 2 + I/6X3, [X 2 , X 6 ] = X 6 , [X 3 , X 6 ] = 0. (65) 

In the three previous subcases, the algebra is ^-dimensional. A 5-dimensional algebra is 
generated in the two remaining subcases. In the subcase n 7^ 2, n 7^ 3 ; a = — 6((4 — 
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2n)/(n - 3) + l), p = -A, the non-rigid dilation X 4 and the temporal inversion X 6 enlarge 
the commutator structure ( I6ip - (l62p fry 

[X 1; X 4 ] = [X 4 ,X 6 ] = 0, [X 2 ,X 4 ] = (2 - n)X 4 , [X 3 ,X 4 ] = 2(n - 2)frX 4 . (66) 

In the subcase n = 2, a = —b, p = —4, the logarithmic dilation X 5 and the temporal inversion 
X 6 enlarge the commutator structure (1611) - (I6"2"j) fry 

[X l5 X 5 ] = [X 5 , X 6 ] = 0, [X 2 , X 5 ] = X 3 , [X 3 , X 5 ] = -2bX 3 . (67) 

(ii) For c ^ 0, p = 1, the point symmetries additionally comprise the shift-scaling X 3 ; the 
non-rigid dilation X 4 when a = — b((l — n/2) /(n — 3) + n 7^ 2, and the logarithmic 

dilation X5 when a = —b, n = 2. In the subcase n ^ 3, d ^ — fr((l — n/2)/{n — 3) + l) ; the 

shift- scaling X 3 enlarges the commutator structure to generate a 3-dimensional algebra fl6Tl) - 
fl62|) . v4 ^-dimensional algebra is generated in the two remaining subcases. In the subcase 
n 7^ 2, n 7^ 3, a = — b((l — n/2)/{n — 3) + l) ; the non-rigid dilation X 4 enlarges the 
commutator structure (!6T!) - (|62|) fry 

[X 1; X 4 ] = 0, [X 2 ,X 4 ] = (2 - n)X 4 , [X 3 ,X 4 ] = (1 - n/2)frX 4 . (68) 

In the subcase n = 2, a = —b, the logarithmic dilation X5 enlarges the commutator structure 
by 

[X 1; X 5 ] = 0, [X 2 , X 5 ] = X 3 , [X 3 , Xb] = fr/2X 3 . (69) 

Variational point symmetries of a radial wave equation (!5a| comprise a sub-algebra of the 
algebra described in Theorem [lj The structure of this sub-algebra will be discussed in the 
context of conservation laws in section 14. 11 

Theorem 2. For a multi-dimensional, hyperbolic, semilinear radial wave equation fl5a|) . its 
point symmetries in the general case n ^ 1, bp — 0, a ^ 0, c + b ^ again comprise a 
2-dimensional algebra generated by the time-translation Xi and the space-time dilation X 2; 
with the commutator structure 

[X 1 ,X 2 ]=X 1 . (70) 

Its additional point symmetries generate larger algebras in following cases: 
(i) For p = 0, b 7^ ; the point symmetries additionally comprise the scaling X 3 ; the non- 
rigid dilation X 4 when (a + b + c)(n — 3) = 0, 11/ 2; and t/ie logarithmic dilation X5 u>/ien 
a + fr + c = 0, n = 2. In the subcase n 7^ 3 ; a 7^ — (fr + c), the scaling X 3 enlarges the 
commutator structure (170!) fry 

[X 1 ,X 3 ] = [X 2 ,X 3 ] = (71) 

generating a 3-dimensional algebra. In the subcases n = 3, or n 7^ 2, n 7^ 3, a = —(b + c), 
the scaling X 3 and the non-rigid dilation X 4 enlarge the commutator structure (I70p and (1711) 
by 

[X 1; X 4 ] = [X 3 , X 4 ] = 0, [X 2 , X 4 ] = (2 - n)X 4 (72) 

generating a ^-dimensional algebra. In the subcase n = 2, a = — (b + c), the scaling X 3 and 
the logarithmic dilation X5 enlarge the commutator structure (!70|) and (!TT|) fry 

[X 1; Xb] = [X 3 , X 5 ] = 0, [X 2 , Xb] = X 3 (73) 
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generating another 4- dimensional algebra. 

(ii) For b = 0, c 0, the point symmetries additionally comprise the non-rigid dilation X 4 
when (a + (1 — p)c){n — 3) = 0, p(p - 1) = 0, n ^ 2; and the shifts X 7 and X 8 when p — 1. 
In £/ie subcases p = 0, n = 3, or p = 0, n 7^ 3 ; a = — c ; i/ie non-rigid dilation X 4 enlarges 
the commutator structure (1701) fry 

[X 1 ,X 4 ]=0,[X 2 ,X 4 ] = (2-n)X 4 (74) 

generating a 3-dimensional algebra. In the subcase p=l, n^3, £/ie s/w/fe X 7 and X 8 enlarge 
the commutator structure (1701) fry 

[Xi, X 7 ] = [X 2 , X 7 ] = 0, [X 1; X 8 ] = X 7 , [X 2 , X 8 ] = X 8 (75) 

generating a ^-dimensional algebra. In the subcase p = 1, n = 3, the non-rigid dilation X 4 
and the shifts X 7 and X 8 enlarge the commutator structure (1701) and (1751) by 

[X l5 X 4 ] = [X 4 , X 7 ] = [X 4 , X 8 ] = 0, [X 2 , X 4 ] = -X 4 (76) 

generating a 5- dimensional algebra. 

Since semilinearity for a radial wave equation (l5al) precludes the existence of a variational 
formulation, the algebra described in Theorem [2] has no variational sub-algebra. 

3.2. Point symmetries of u tt = (c + (a + b)u p r )u rT + (n — l)(cu r + bu p r +l )/r. 

Table [2] lists the infinitesimal point symmetries of the radial wave equation (15bj) in all cases 
such that this equation is nonlinear (i.e. pb ^ or p(a + b) 7^ 0), hyperbolic (i.e. p{b + a) 7^ 
or c 7^ 0), and multi-dimensional (i.e. n 7^ 1, c 7^ or (p + 1)6 7^ 0). These three restrictions 
are equivalent to the inequality (n — l)p(b 2 + c 2 )((6 + a) 2 + b 2 c 2 (p + l) 2 ) 7^ holding on the 
nonlinearity power p, the dimension n, and the coefficients a, b, c. 

Symmetry X 4 generates a time-translation. Symmetry X 2 generates a scaling on (t,r,u). 
Symmetries X3 and X 4 generate a shift on u. These are the only symmetries admitted for 
all allowed values of a, 6, c, p, n. Symmetry X 5 generates a scaling on (t, u), while symmetry 
X6 generates a temporal inversion on (t,u). Symmetry X 7 is a non-rigid generalization of 
the scaling on r given by X 2 + |Xs (i.e. r = 0, £ = r, r\ = 0) when c = 0. Symmetry X 8 is a 
logarithmic counterpart of X 7 . 

From table [2] we see that the point symmetry structure for the wave equation (15bl) is 
richest in the case c = 0, i.e. when the linear Laplacian term is absent, and in the subcases 
p = —2 and p = —4, i.e. when the nonlinearities take the form 

Utt = ^( U ;\% + (n-l)bu;\ p+1 \/ r , p = -2,4. (77) 

IpI 

We also see that there is no distinguished dimension n / 1, while the only distinguished 
relation among a, b, c is given by 

a = b(n-2), b^0. (78) 

When the radial wave equation ( |5bl) is quasilinear (i.e. p 7^ 0, b + a 7^ 0), this relation (17H1) 
coincides with the condition a = pb for existence of both a variational structure and a 
divergence structure ([8]) ifp = n — 2^0. Moreover, the relation (178|) corresponds to 

m = p + 1 (79) 
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in terms of the radial power 



which characterizes the alternative variational and divergence structures (T16l) and ( 1151) that 
exist for the radial wave equation (l5bl) in the highly quasilinear case c = (p 7^ 0, 6 + a 7^ 0). 
In contrast, when the radial wave equation fl5b|) is semilinear (i.e. p 7^ 0, b + a = 0, he 7^ 0), 
the non-existence of both a variational structure and a divergence structure precludes any 
overlap with the relation ([75]) . 

We will next summarize the structure of the algebra generated by all of these point sym- 
metries. 



Theorem 3. For a multi-dimensional, hyperbolic, quasilinear radial wave equation ( I5bj) . its 
point symmetries in the general case n / 1, 6 / 0, 6 + fl / 0, p(p + 1) 7^ comprise a 
4- dimensional algebra generated by the time-translation X 1; the scaling X 2 , and the shifts X 3 
and X 4; with the commutator structure 

[X 1; X 2 ] = Xx, [X 1; X 3 ] = [X 2 , X 4 ] = [X 3 , X 4 ] = 0, [X 2 , X 3 ] = -X 3 . (81) 

For c = 0, its point symmetries generate larger algebras that additionally comprise the tem- 
poral scaling X 5 ; the temporal inversion X 6 when p = —4; the non-rigid scaling X 7 when 
a 7^ b(n — 2), p = —2; and the logarithmic scaling X 8 when a = b(n — 2), p = —2. In the 
subcase p 7^ —2, p 7^ —A, the temporal scaling X 5 enlarges the commutator structure ( IHTj) by 

[X 1; X 5 ] = pX 1; [X 2 , X 5 ] = 0, [X 3 , X 5 ] = -2X3, [X 4 , X 5 ] = -(2 + p)X 4 (82) 
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generating a 5- dimensional algebra. In the subcase b ^ 0, p = —A, the temporal inversion X 6 
enlarges the commutator structure (IHTj) and (1821) by 

[Xi,X 6 ] = — |X 5 , [X 2 ,X 6 ] = X 6 , [X 3 ,X 6 ] = X 4 , [X 4 ,X 6 ] = 0, [X 5 ,X 6 ] = — 4X 6 . (83) 

In the subcase a ^ b{n — 2), p = —2, the non-rigid scaling X 7 enlarges the commutator 
structure flHTT) and fl82l) by 

[X X ,X 7 ] = [X 3 ,X 7 ] = [X 4 ,X 7 ] = [X 5 ,X 7 ] = 0, [X 2 ,X 7 ] = -(m + 1)X 7 . (84) 

In the subcase a = b{n — 2), p = —2, the logarithmic scaling X 8 enlarges the commutator 
structure (I8T1) and (1821) by 

[X 1; X 8 ] = [X 3 , X 8 ] = [X 4 , X 8 ] = [X 5 , X 8 ] = 0, [X 2 , X 8 ] = X 2 + ±X 5 . (85) 
In the previous three subcases, the algebra is 6- dimensional. 

Variational point symmetries of a radial wave equation f!5bj) comprise a sub-algebra of the 
algebra described in Theorem [3J The structure of this sub-algebra will be discussed in the 
context of conservation laws in section 14.21 

Theorem 4. For a multi-dimensional, hyperbolic, semilinear radial wave equation (15bl) . its 
point symmetries in all cases n ^ 1, 6 + fl = 0, 6 / 0, c ^ 0, p(p + 1) ^ comprise a 
^-dimensional algebra generated by the time-translation X 1; the scaling X 2 , and the shifts X 3 
and X 4 , with the commutator structure 

[Xi,X 2 ] = Xi, [X 1; X 3 ] = [X 2 ,X 4 ] = [X 3 ,X 4 ] = 0, [X X ,X 3 ] = X 3 , [X 2 ,X 3 ] = -X 3 . (86) 

Since the existence of a variational formulation is precluded by semilinearity for a radial 
wave equation (|5b]), the algebra described in Theorem [4] has no variational sub-algebra. 

3.3. Computational Remarks. 

In the computation of point symmetries (1191) for a radial wave equation ([5]), the symmetry 
determining equation f|2Up is formulated in the jet space whose coordinates are defined by 
t, r, u, u t , u r , u tr and u rr , where u u is replaced by gu„ + f. For wave equation fl5aj) . this 
determining equation splits with respect to u t , u r , u tr , u rr into an overdetermined system 
of 16 equations in 8 unknowns r](t,r,u), £(t,r,u), r{t,r,u), a, b, c, n, and p. In contrast, 
for wave equation f )5b|) . the determining equation splits initially with respect to u t , u tr and 
u rr , but not with respect to u r because this wave equation contains terms involving the 
unknown power p of u T . In this case the partial splitting yields an overdetermined system of 
9 equations in the 8 unknowns r](t, r, u), £(t, r, u), r{t, r, u), a, b, c, n, and p. 

We use computer algebra both to carry out the splitting and to solve the resulting overde- 
termined system of equations. The formulation and splitting of the symmetry determining 
equation is done automatically by the program LiePDE [17J which in turn calls the package 
Crack [19J to solve the resulting overdetermined system for the 8 unknowns. Although 
the system is linear in the unknowns rj, t and £, it is nonlinear jointly in the 8 unknowns 
and even is non-polynomial in the unknown power p which appears as an exponent. Such 
problems are typically much more complicated to solve than are linear problems. 

One main difficulty is the handling of unknowns as exponents, which requires case splittings 
such that all possible balances of exponents are considered. An example is the equation 
£ t — (a + b)uP.r r — cr r — (a + b)u^. +1 r u — cu r r u = together with the inequality p ^ arising 
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during the computation for solving the overdetermined system in the case of wave equation 
(|5b|) . Balancing the powers 0, 1, p, p + 1 of u r leads to three different splittings: p = —1; 
p = 1; (p — l)(p + 1) 7^ 0. Each of these three cases then must be solved separately in 
the remaining computation. A similar situation arises during the computation for solving 
the overdetermined system for wave equation (IBli|) . which involves splitting with respect to 
powers of u after integrating out the u dependence of r, £, r\. 

Other kinds of case splittings also arise. For example, when an expression factorizes then 
each factor must be considered zero or non-zero, as in the equation c(£ r — r t ) = arising 
during the computation for solving the overdetermined system for wave equation f!5bl) . This 
leads to the cases c = and c ^ 0, £ r = Tj, both of which must then be solved separately. 

To our knowledge, Crack is the only package that can handle such computational diffi- 
culties and, in particular, can both generate and solve all resulting cases automatically. In 
handling nested case distinctions, Crack considers equations and inequalities together. For 
example, equations and inequalities are simplified with each other, and new inequalities as 
well as equations are generated through case distinctions. 

For both wave equations (jSaj) and f!5bl) . Crack is able to complete the entire computation 
of point symmetries automatically. 



For each class of radial wave equations ( |5a|) and f|5bl) we will now find all kinematic conser- 
vation laws fl56l) - fl57j) and all generalized-energy conservation laws (j4"6]) - (j4"T)) . In particular 
we will explicitly determine any such conservation laws that exist only for special nonlin- 
earity powers p and dimensions n ^ 1, as well as for special relations among the constant 
coefficients a, b, c in these wave equations, excluding cases where the wave equation is linear 
or non-hyperbolic. (Note we will allow n to have non-integer values. An interpretation of 
the wave equations (l5a|) and (|5b|) in such cases is given in section [5j) 

Our results are obtained by first solving the determining system (|30|) - (!3T1) for the multiplier 
function 



given by the characteristic form f !2B|) of these conservation laws, and then using the homo- 
topy integration formula ( l3~4"|) - f[3T)j) or directly integrating relations f l3"!?j) -( l3l?j) to obtain the 
conserved density T(t,r,u,u t ,u r ) and flux X(t,r,u,u t ,u r ) determined by each multiplier. 
To keep the list of solutions succinct, we have merged conservation laws with similar forms 
into a combined form wherever possible. Remarks on the computation will be provided at 
the end in section 14.31 

In the case when a radial wave equation (I5al) or (15b j) possesses a variational formulation 
(|9]), multipliers (!H7|) correspond (via Noether's theorem (!39l) - (H0|) ) to symmetries with the 
characteristic form 



Thus, our classification of all generalized-energy conservation laws will provide a partial 
classification of variational symmetries for these wave equations. We note that a symmetry 
generator ( 188]) will be of point type ( !T9|) iff 



4. Conservation law Classification 







da d 2 /3 
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du r du r du r 
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Table 3. Kinematic conservation laws for 1, (bp) 2 + a 2 (c + 6) 2 7^ 0. 



in which case it will have the equivalent canonical form 



X = -r 1 



ad/di + --d/dr + (^-u r - (3)0/ du) . (90) 



In the case of the alternative variational structure f lT6|) - f|T71) . the same statements hold with 

r l-n re pl acec l py r ~m_ 



4.1. Conservation laws of u u = (c + bu p )(u rr + (n — l)u r /r) + au p ~ x u 2 . 

In the following tables we list the kinematic conservation laws (l56l)- (|571) and generalized- 
energy conservation laws (H6l) -()471) of the radial wave equation (l5al) in all cases such that 
this equation is nonlinear (i.e. a 7^ or pb 7^ 0), hyperbolic (i.e. pb 7^ or c + 6 / when 
pb = 0), and multi-dimensional (i.e. n / 1). These three restrictions are equivalent to the 
inequality (n — l)((fyo) 2 + a 2 (c + b) 2 ) 7^ holding on the nonlinearity power p, the dimension 
n, and the coefficients a, b, c. 

To begin, the kinematic conserved densities, fluxes, and corresponding multipliers are listed 
in Table El where H(u) denotes the expression ffT31) appearing in the divergence structure 

(TO. 

Conserved densities T± and T2 arise from the n-dimensional divergence structure (jHJ) for the 
radial wave equation (|5ap . Conserved densities T 3 and T 4 arise from an alternative divergence 
structure {ru t )t = ijF) r where the flux is given by F = (c + bu p )u r + (n — 2) (cm + bH{u)), 
which holds under the same existence condition a = pb as for the divergence structure flH]). 
A logarithmic counterpart of the alternative divergence structure gives rise to the conserved 
densities T 5 and T 6 . 

We will now summarize the structure of the vector space spanned by the conserved quan- 
tities ( 12 5 p arising from these kinematic conserved densities. 

Theorem 5. A multi- dimensional, hyperbolic, nonlinear radial wave equation (15 ap admits 
kinematic conservation laws only in cases when a = pb 7^ 0. 
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Table 4. Generalized-energy conservation laws for 1, (6p) 2 + a 2 (c + 6) 2 ^ 0. 



(i) For n ^ 2, a = p6 ^ 0, a basis of kinematic conserved quantities is comprised by 

poo poo poo poo 

Ci= u t r n ~ l dr, C 2 = (tu t - u)r n ~~ l dr, C 3 = u t rdr, C 4 = {tu t - u)rdr. 
Jo Jo Jo Jo 

(91) 

(ii) For n = 2, a = pb ^ 0, a basis of kinematic conserved quantities is comprised by 

poo poo poo poo 

C\ = / u t rdr, C 2 = (tu t —u)rdr, C 5 = / u t (\nr)rdr, C 6 = / (tu t —u) (In r)rdr. 
Jo Jo Jo Jo 

(92) 

In both cases the conserved quantities span a ^-dimensional vector space. 

The conditions on the coefficients a, 6, c coming from Theorem [5] combined with Proposi- 
tion [1] imply the following interesting result. 

Corollary 1. A multi-dimensional, hyperbolic, nonlinear radial wave equation fl5a|) possesses 
kinematic conservation laws iff it has a quasilinear divergence structure. 

Next, the generalized-energy conserved densities and fluxes are listed in Table HI corre- 
sponding multipliers are listed in Table |5j 

From Table |4] we see that these conservation laws are characterized by two distinguished 
relations among the coefficients a, b, c: 2a = pb, or 3a = pb. The first relation 2a = pb is 
simply the condition for existence of a variational structure ([9]) for the radial wave equation 
( )5al) . Conservation laws in this case correspond to variational symmetries through Noether's 
theorem, as we will summarize in the next subsection 14.1.11 The second relation 3a = pb turns 
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Q 


conditions 


1 


r Ut 


2a = 6p 


2 


r" _i (b(np + 4)tu t + b(p + A)ru r 
+ 2{n — l)(c + on) J 


2a = 6p, c(p — 1) = 


3 


r(bru r + bu + c) 


2a = bp, c(p — 1) = 0, 
p = 2(n - 3) 


4 


r n -H{tu t - u) 


2a = ftp, c = 0, 
p = —4 


5 


rtut + (r In r){ru r + u) 


2a = 6p, c = 0, 
p = — z, n = z 


6 


r n ~ 1 u t (36ru r + + c)) 


3a = ftp, c(p — 1) = 0, 
p = 6(1 — 3/n) 


7 


r 4 ' 5 (ttt t — u)(5ru r + 3u) 


3a = 6p, c = 0, 
p = —4, n = 9/5 



Table 5. Generalized-energy multipliers for n ^ 1, (fyo) 2 + a 2 (c + 6) 2 ^ 0. 



out to be the condition needed for deriving a variational type of Morawetz radial dilation 
identity which produces conservation laws from temporal symmetries. We will discuss this 
interesting case later in subsection 14.1.21 

4.1.1. Noether correspondence for u tt = (c + bu p )(u rr + (n — l)u r /r) + (bp/2)u p ~ l u1. In the 
variational case 2a = pb, multipliers Q for conservation laws of the radial wave equation 
(|5ap are in one-to-one correspondence to variational symmetries X = r x ~ n Qdjdu of the 
Lagrangian fllOap . We see from Table that all of the generalized-energy multipliers for 
2a = pb are linear in u t and u r . Consequently, the corresponding variational symmetries 
are all of point type. This correspondence is summarized in Table El where C denotes 
the conserved quantity (125]) given by the conservation laws, and X denotes the variational 
symmetry generator in canonical form fl58|) . 

Conserved quantity C\ is an energy arising from the time-translation symmetry X x . Con- 
served quantity C<i is a similarity energy arising from the scaling symmetry X 2 , while con- 
served quantity C5 is a logarithmic counterpart of C2. Conserved quantity C3 is a radial 
momentum arising from the non-rigid radial dilation symmetry X3. Conserved quantity C4 
is a conformal energy arising from the temporal inversion symmetry X4. 

We will now state the vector space structure spanned by these conserved quantities as well 
as the algebra structure determined by the corresponding variational symmetries. 

Theorem 6. For a multi-dimensional, hyperbolic, nonlinear radial wave equation f!5a|) that 
has a variational structure (Q when 2a = pb ^ 0, variational point symmetries are admitted 
in two cases, corresponding to all generalized- energy conservation laws. 
(i) Fore = 0, the variational point symmetries comprise the time-translation Xi/ the scaling 

X 2 ; the non-rigid radial dilations X 3 when p = 2(n — 3) and X 5 when p = —2, n = 2; and 
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X 


C 


conditions 


1 


-d/dt 


POO 

/ l(u 2 t + (c + bu p )u 2 r )r n - 1 dr 
Jo 




2 


— h ( ti T) -\- A\tf) f)t — h( 11 -\- A\rr) I r)r 

U\ILU \^ rt J l \J / \J b "\ A' ) {J f {J > 

+ 2(n- l)(c + bu)d/du 


I'iXJ 

/ \ (b(np + A)t(u 2 t + (c + bu p )u 2 r ) 
Jo 

+ 2b{p + 4)ru r u t 

+ 4(n - l)(c + 6M)ut)r n_1 dr 


c(p- 1) = 


3 


—br 3 ~ n d / dr + r 2 ~ n {bu + c)d/du 


/•oo 

/ (&rw r + 6m + c)w t r dr 
io 


c lP — X J — U ; 
p = 2(n - 3) 


4 


—to/ at — tuo/ou 


/'OO 

/ l(t 2 (u 2 + bu- 4 u 2 ) 
Jo 

- 2tuu t + u 2 )r n ~ l dr 


c = 0, 
p = -4 


5 


—td/dt + r In rd/dr + (lnr)ud/du 


/'0<J 

/ (\t{u 2 t + bu- 2 u 2 .) 
Jo 

+ (In r)(ru r + u)u/jr dr 


c = 0, 

p = -2, n = 2 



Table 6. Generalized-energy conserved quantities and variational symmetries 
for n ^ 1, 2a = pb ^ 0. 



£/ie temporal inversion X 4 u>/ien p = —4. In t/ie subcase p 7^ 2(n — 3), p 7^ — 4, t/ie tzme- 
translation Xi and t/te scaling X 2 generate a two-dimensional algebra with the commutator 
structure 

[X 1 ,X 2 ] = -6(np + 4)X 1 . (93) 

TTie corresponding two-dimensional vector space of conserved quantities is spanned by the 
energy C\ and the similarity energy C^- In the subcase p = 2(n-3), p ^ —4, n 7^ 2, the 
non-rigid radial dilation X 3 enlarges the commutator structure ( 193|) by 

[X 1; X 3 ] = 0, [X 2 ,X 3 ] = 2b(n - l)(n - 2)X 3 (94) 

generating a three-dimensional algebra, while the corresponding vector space of conserved 
quantities is spanned by C\, C 2 , and the radial momentum C 3 . In the subcase p = —A, the 
temporal inversion X4 enlarges the commutator structure ( 193]) by 

[X 1; X 4 ] = -1/(46(72 - 1))X 2 , [X 2 ,X 4 ] = Ab(n - 1)X 4 (95) 

generating another three-dimensional algebra, while the corresponding vector space of con- 
served quantities is spanned by C\, C 2; and the conformal energy C4. In the subcase p = —2, 
n = 2, the non-rigid radial dilations X 3 and X5 enlarge the commutator structure (|93|) - (|94l) 
by 

[X 1} X 5 ] = -X x , [X 2 , X 5 ] = 2X 3 , [X 3 , X 5 ] = X 3 (96) 

generating a four-dimensional algebra whose corresponding vector space of conserved quan- 
tities is spanned by C\, C 2 , C 3; and the dilational energy C5. 

(ii) For c 7^ 0, p = 1, the variational point symmetries comprise the time-translation Xi; the 

scaling X 2 ; and the non-rigid radial dilation X 3 when n = 7/2. In the subcase n 7^ 7/2, the 
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time-translation and the scaling X 2 generate a two-dimensional algebra with the commu- 
tator structure (1931) . while the corresponding vector space of conserved quantities is spanned 
by the energy C\ and the similarity energy CV In the subcase n = 7/2, the non-rigid radial 
dilation X 3 enlarges the commutator structure f[9"31 by 

[X 1 ,X 3 ]=0,[X 2 ,X 3 ] = 15/26X 3 (97) 

generating a three-dimensional algebra, while the corresponding vector space of conserved 
quantities is spanned by C\, C2, and the radial momentum C3. 

4.1.2. Morawetz-type variational identity for u u = (c+bu p )(u rr + (n — l)u r /r) + {bp/?>)u p ~ l u1. 
To derive a variational identity for the radial wave equation (I5a|) when 2a ^ pb, we start by 
making use of the previous Lagrangian fllOap to write 

BL 

r n ^(u tt - (c + bu p )(u„ + (n- l)u r /r) - au p - l u 2 r ) = E u (r n ~ l L) + /ir"" 1 — (98) 

Ou 

where 

2a , , . 

/i = l--^0 (99) 

is the coefficient of the non-variational terms, and where 

E u = 3/du - D r d/du r - Dtd/dut (100) 

is the 1st order Euler-Lagrange operator (i.e. local variational derivative) with respect to 
u. We next multiply equation ( 198]) by the expression W = ru r — Am which represents the 
characteristic function given by the generator of a scaling transformation on (r,u), namely 
St = 0, Sr = r, Su = Xu, for some constant A. This yields 

r n ~ 1 (rw r - Am) (u tt - (c + bu p )(u rr + (n - l)u r /r) - avP^ul) 

dL 

= (ru r - \u)EJr n - l L) + ^r n - l (ru r - Am)—. (101) 

ou 

Now we will show that the Euler-Lagrange term in this equation can be expressed in a 
useful way as the sum of a total variational derivative i? M (r n_1 L) for some new Lagrangian 
±j[U } M r , V>i ) and a dilational derivative (rD r + A)(r n 1 f) for some function f(u, u t ) and some 
constant A. Our main tool is the operator identity 

WE U = prY - D t (Wd/du t ) - D r (Wd/du r ) (102) 

where prY is the prolongation of the generator Y = Wd/du. 

To proceed we split up the Lagrangian L = Lq + L\ into its temporal and radial parts 

L = -\u\, Li = |(c + bu p )u 2 r (103) 

which have the scaling properties 

prYL = rD r L - 2AL (104) 

and 

dL 

prYL x = rD r L 1 - 2(A - l)L x - Am— - (105) 

ou 

for W = ru r — Am. Applying the identity (I102j) to the part of the Euler-Lagrange term 
containing L in equation fllOip . and using the scaling property (11041) . we get 

(ru r - \u)E u {r n - l L ) = A(r n-1 (rM r - Xu)u t ) + r n ~\rD r - 2X)L . (106) 
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If we now define Lq by 

dL 



(ru r — Xu)u t (107) 



du t 

then the time derivative term in equation (I106p can be expressed as 

D t {r n " l {ru r - \u)u t ) = -S tt (r n_1 L ) + r n ~\rD r + (A + n))L (108) 

where 

L = \(ru r - \u)u 2 t . (109) 

Hence equation (11061) becomes 

(ru r - X^E^r^Lo) = -E u {r n ~ l L ) + r n -\2rD r + (n - A))L . (110) 

Similarly, from the scaling property (11051) combined with the identity (I102p applied to the 
part of the Euler-Lagrange term containing L\ in equation (11011) . we obtain 

(ru r — \u)E u (r n ~ 1 Li) = -D r (r n -\ru r -\u)(c+bu p )u r )+r n - l ((rD r -2(\-l))L 1 -\u—±). 

ou 
(111) 

If we integrate by parts on the dilational derivative term in equation fillip and then define 
Li by 

dL 

' -\{ru r -2\u){c + bu p )u Tl (112) 



du 



r 



we find that the radial derivative terms in equation fillip can be expressed as 

BL 

-D r (r n -\ru r -\u)(c+bu p )u T )+r n D r L 1 = -E M (r n - 1 L 1 )+r n - 1 ((A(p+l)-n)L 1 -|ra r ^) 

(113) 

where 

l x = -l(ru r -3\u)(c + bu p )u 2 r . (114) 
Hence equation (II 111) becomes 

(rUr-X^E^ 1 ^) = -E u (r n - 1 L 1 )+r n - 1 ((X(p-l) + 2-n)L 1 -(lru r + Xu)-^). (115) 

Finally, substituting expressions (IllOp and f 1 1 1 5 j) into equation (UOip . we have 
r n_1 (ru r — Xu)(u tt — (c + bu p ){u Tr + {n — l)u r /r) — au^u 2 ,) 
= -E u (r n -\L + Lt)) + r n -\2rD r + (n - X))L 

81 

+ r 71 - 1 ((A(p - 1) + 2 - n)L x + ((/i - \)ru r - (1 + ^)Au)— ^) • (116) 

We now require that the terms Li, u r dLi/du, udL\/du vanish, which holds iff the coefficients 
of u p ~ l ul, u 2 ., u p u 2 in this equation are equal to zero. This imposes the respective conditions 

^ = 1/3, cA = 0, (l+p/3)A = 2-n. (117) 

From the expression f )99|) for fi in terms of a, b, we then obtain 

3a = pb. (118) 
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Thus, for 3a = pb, we see that equation (jlOip reduces to the variational dilation identity 
r n ~ x {ru r - Xu) (u u - (c + bu p )(u rr + (n - l)u r /r) - (p6/3)M p_1 ^) 
= -Euir^L) - {rD r + A)(r n "\ 2 ) (119) 



with 
and 



L + L 1 = \{ru r - Xu)u\ - \{ru r - 3Xu)(c + bu p )u 2 r (120) 



A = l-(n + A)/2. (121) 
As we will now illustrate, this identity f ll 19j) is able to produce conservation laws from 
multipliers that are given by certain symmetries X = Pd/du in characteristic form, in a 
manner similar to Noether's theorem (l39l)-(H0i). 
Consider the time translation 

X = u t d/du. (122) 
When the identity ( 1119ft is multiplied by 

P = ut (123) 

we see first that the Euler-Lagrange term yields a total time derivative 

8T ~ BT 

- u t E u (r n ~ l L) = D t (r n - l u t —) - r^prXL = D t (r n - l u t — - r^L) (124) 

since, by inspection, the Lagrangian functional C = J °° Lr n ~ 1 drdt is invariant under 
time-translation. Second, we see that the dilational derivative term produces 

-utirDr + X)^' 1 ^) = -2r n u 2 t u tT -{X + n-l)r n ~ l u] = D r (- \r n ul) + (1 - \n- X)r n - l u 3 t 

(125) 

which will reduce to a total radial derivative iff 

X = l-n/3. (126) 

Equating expressions A126D and (1121ft . and using conditions ( 1117ft . we obtain the additional 
conditions 

c = 0, p = 6(l-3/n) (127) 

as well as 

A = -n/3. (128) 

Thus, for a, b, c, p, X satisfying conditions ( 1118ft . ( 1126ft and ( 1127ft . the product of the 
variational dilation identity ( 1119ft and the symmetry characteristic function ( 1123ft yields a 
conservation law, where the conserved density and flux are given by T 6 and X 6 in Table H] 
while the multiplier is given by Qq in Table El We note that the form of this multiplier 
consists of r n ~ l times the product of the characteristic functions W = ru r — Xu and P = u t 
representing the scaling transformation and the time-translation symmetry respectively. The 
c / case for this conservation law shown in Tables H] and [5] can be obtained from a 
similar variational identity in which the scaling transformation is generalized to a shift- 
scaling transformation given by a characteristic function of the form W = ru r — X(u + v), for 
some constants A 7^ and v 7^ 0. In addition, the remaining conservation law given by the 
conserved density T 7 and flux X7 in Table HI with the multiplier given by Q7 in Table can 
be obtained from the variational identity (1119ft by using the temporal scaling symmetry with 
the characteristic function P = tu t — u in place of the time-translation symmetry. These 
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Y 


X 


C 


conditions 


6 


n(bu + c)d/du 
— 3brd/dr 


—a/at 


/ | ((36rw r + n(bu + c))w 2 

+ (bu p + c)(bru r + n{bu + c))w 2 )r n_1 dr 


c(p-l) =0, 
p = 6(1 — 3/n) 


7 


3ud/du — 5rd/dr 


—ud/du — td/dt 


/'(X) 

/ | (u r (5bru r + Qbu)(u — tu t )u~ 4 
Jo 

+ \ru 2 t {3u -iu t ))r 4/5 dr 


c = 0, 

p = —4, n = 9/5 



Table 7. Generalized-energy conserved quantities and generators for n^l, 
3a = bp 7^ 0. 



results are summarized in Table where C denotes the conserved quantity (1251) given by 
the conservation laws, and where Y and X respectively denote the scaling generator and 
symmetry generator in canonical form ( 1581) . 

We will now state the vector space structure spanned by these conserved quantities. 

Theorem 7. For a multi-dimensional, hyperbolic, nonlinear radial wave equation (|5a|) that 
has no variational structure (19~1), generalized- energy conservation laws are admitted only in 
cases when 3a = pb ^ 0, arising from a variational type of radial dilation identity. 

(i) For c = 0, 3a = pb 0, the conserved quantities are comprised by the generalized- 
dilation energy Cq when p = 6(1 — 3/n) and the generalized- dilation conformal energy C-j 
when p = —4, n = 9/5. In the subcase p = 6(1 — 3/n), n ^ 9/5, a basis for the conserved 
quantities consists of Cq, spanning a 1-dimensional vector space. In the subcase p = —4, 
n = 9/5, the basis consists of Cq, C?, spanning a 2-dimensional vector space. 

(ii) For c ^ 0, 3a = pb ^ 0, p = 1, the conserved quantities comprise only the generalized- 
dilation energy Cq when n = 18/5, spanning a 1-dimensional vector space. 

4.2. Conservation laws of u tt = (c + (a + b)u^)u rr + (n — l)(cu r + bu p+l )/r. 

In the following tables we list the kinematic conservation laws (15T)|) - (15"T1) and generalized- 
energy conservation laws (r46l) - (T4TI) of the radial wave equation (I5bj) in all cases such that 
this equation is nonlinear (i.e. p{p + 1)6 ^ or p(a + 6) ^ 0), hyperbolic (i.e. p(6 + a) ^ 
or c ^ 0), and mult i- dimensional (i.e. n ^ 1, c ^ or (p+ l)b ^ 0). These three restrictions 
are equivalent to the inequality {n — l)p(6 2 + c 2 )((6 + a) 2 + 6 2 c 2 (p + l) 2 ) ^ holding on the 
nonlinearity power p, the dimension n, and the coefficients a, 6, c. Wherever it is convenient, 
we will use the notation 

k = a + b (129) 

and 

m = b(n- l)(p + l)/(o + 6). (130) 

To begin, the kinematic conserved densities, fluxes, and corresponding multipliers are 
listed in Table [H 

These conservation laws are characterized by three distinguished relations among 6, c, k: 
k = (p + 1)6 (i.e. a = pb), or c = 0, or k = (p + l)(n — 1)6 (i.e. a = {p{n — 1) + n — 2)6), 
c ^ 0. When k = {p + 1)6 the radial power (I130p reduces to m = n — 1, so thus in this 
case the conserved densities Ti and T 2 arise from the n-dimensional divergence structure ([H]) 
for the radial wave equation (15bj) with c ^ 0. When c = the conserved densities T\ and 
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T 


X 


conditions 


1 






r m (cu r + b(yU 1 Kp +1 ) 
\ m J 




r(k — (n+ l)b) = 


2 


r m t 


r m (tu t — u) 


r m t(cu r + b(yT} '~ ' X K p r +l ^ 
\ m / 


) 


c(k- (p + l)b) = 


3 


r 


ru t 


—b(n — l)ru p+1 — c[ru r + in — 


2)u) 


k= (p + l)(n-l)6,c^0 
n ^ 2 


4 


rt 


r(tut — u) 


—t(b(n — l)ruP +1 + c(ru r + (n - 


-2)u)) 


k = {p + l){n - 1)6, c / 
n ^ 2 



Table 8. Kinematic conservation laws for [n — l)(6 2 (p + l) 2 + c 2 ) 7^ 0, p(k 2 + 



6 2 c 2 (p+l) 2 ) 7^0. 



T2 arise from the alternative divergence structure ( Tl5i) for this wave equation ( I5b[) . with an 
arbitrary radial power (I130p . In contrast, when c^O and k = (p+ l)(n — 1)6, the conserved 
densities T3 and T4 arise from a different divergence structure (rut)t = {i r F) r where the flux 
is given by F = c(ru r + (n — 2)tt) + (k/(p + 1))m^ +1 . We note that, in the case c = 0, 
this divergence structure coincides with the alternative divergence structure (|T5|) . since the 
relation k — (p + l)(n — 1)6 is equivalent tom=l. 

We will now summarize the structure of the vector space spanned by the conserved quan- 
tities ( 12 5 p arising from these kinematic conserved densities. 



Theorem 8. A multi-dimensional, hyperbolic, nonlinear radial wave equation (15bj) admits 
kinematic conservation laws only in the following cases: 

(i) For c^O, k = (p + 1)6 7^ 0, a basis of kinematic conserved quantities is comprised by 



C x = / u t r n - L dr, C 2 = / (fa ~ u)r n ~ l dr. (131) 
Jo Jo 

(ii) For c 7^ 0, k = (p + l)(n — 1)6 7^ 0, n 7^ 2, a basis of kinematic conserved quantities is 
comprised by 



C3 = u t rdr, C 4 = / (tu t — u)rdr. (132) 
Jo Jo 

(iii) For c = 0, k 7^ (p + 1)6, a basis of kinematic conserved quantities is comprised by 

roo /*oo 

d= u t r m dr, C 2 = (tu t - u)r m dr. (133) 



In all cases the conserved quantities span a 2-dimensional vector space. 

The conditions on the coefficients b,c,k = a + 6 coming from Theorem [8] combined with 
Proposition [Tj imply the following interesting result. 

Corollary 2. A multi-dimensional, hyperbolic, nonlinear radial wave equation ( j5bj) possesses 
kinematic conservation laws iff it has a quasilinear divergence structure. 

Next, the generalized-energy conserved densities and fluxes are listed in Table [9J corre- 
sponding multipliers are listed in Table [TDJ 
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T 


X 


conditions 


1 


r n -\\(u 2 f +cu 2 r ) + bh(u r )) 


—r n ~ x {cu r + bu p r +l )u t 


k = (p + 1)6, c ^ 


2 


r~ 2 (ft(u 2 + cu 2 + 2bh(u r )) 
+ (ru r - u)u t ) 


-\r~'*(r(u A + cu 2 + 2kh(u r )) 

+ 2(c + bu*)(tut - u)u r ) 


k = (p + l)6,c ^ 

77 = —1 


3 


r™U u 2 + 6(n_1) /i(wr)) 


—r m (kuth(u r )) 


C= 


4 


±r m ((p(m + 3) + 4)t(w 2 + — ^— — -h{u r )) 
+ 2np + 4)ru r + {2m — p)uju t ) 


-ir m ((p + 4)r(u 2 + 2/e/i(u r )) 
+ 2((p(m + 3) +4)tw t 

+ (zm PjMl n r i 

777, 


c = (J 
p^ -4 


5 


\r m (3(tu t - u)u t + ktu~ 2 ) 


\r m (k(2tu t - u)u- 3 ) 


c = 0, 
p = -4 


6 


lr m (3(tu t - uf + kt 2 u~ 2 ) 


lr m (kt(tu t -u)u; 3 ) 


c = 
p = -4 


7 


U r Ut 


— \{u 2 + 2Hnw r ) — (n — 1)6 In r 


c = U, 
p = -2 


8 


1 —21 In \\ 

[u t (2ru r — u)) 


-\r~ 2 {b(n - l)uu~ 3 

+ r(2u 2 -36(n- l)u; 2 )) 


c = U, 

fc = 30(77 — 1 )/ 2, 

p = -4 


9 


\r~ x ((2m — tu t )u t + 2(r In r)u r u t 
+ 2H In u r ) 


-|((A;(2(1 + lnw r ) + lnr) + w 2 ) lnr 
+ 2k(tu t - w)«~V _1 ) 


^ — ("1 

c — u, 

A; = 6(77 — 1), 
p = -2 


id 


\2 u t ~ k{i. — mu r ))u r — b[n — L)u/r 


— g + ok in w r Jiit 


c = 0, 
p = -2 


n 


lr m (u(u 2 t + 2b{n ~ 1} h(u r )) 
V m 

— r(2k / h(u r )du r + u r u 2 )) 

V + 2 J / 


r m (^|r(|«? + *M«r)K 
& ( n - 1) P +i\ 

UU t U^ ) 

777 


c = 0, 

4k = 6(1 -n)(p-2), 
p^-2 


12 


|r _2 (2(2u — tu t )(u — 2>ru r )u t 

+ bt(n — l)(9ru r — u)u~ 2 ) 


ir" 2 (6(r7 - l)(9ru r - 2u){tu t - u)u~ 3 
— 2ru 2 {tu t — 3m)) 


c = 0, 

4A; = 6(1 - 77)(p- 2), 
p = -4 



Table 9. Generalized-energy conservation laws for (n— l)(6 2 (p+l) 2 + c 2 ) 7^ 0, 
p(A; 2 + 6 2 c 2 (p + l) 2 ) ^0. 



We see from Table [9] that these conservation laws are characterized by two distinguished 
relations among the coefficients 6, c, k: c = 0, or k = (p + 1)6, c ^ 0. These relations 
are simply the conditions for a radial wave equation (15b j) to possess an n- dimensional vari- 
ational structure OH]) or an alternative variational structure (I16p based on the radial power 
( I128p . Thus, all generalized-energy conservation laws for this wave equation correspond to 
variational symmetries through Noether's theorem. We will summarize this correspondence 
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conditions 
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n — 1 


k = (n + l)b r ^ 
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r~ 2 {tu t + ru r — u) 


k = (p + l)6,c ^ 
n = — 1 
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r m u t 


c = 


4 


r m ((p(m + 3) + A)tut 

-r \P -r QjrUr -+- [zm — p)u j 


c = 
p^ -4 


5 


r m (2tu t - u) 


c = 0, 
p = -4 


6 


r m t(tu t - u) 


c = 
p = -4 


7 


u r 


c — u, 
p = -2 






c = U, 


o 
O 


r~ 2 (2ru r — u) 


p = -4 






— n 


9 


r~ 1 (r(lnr)u r — tut + u) 


A; = b(n — 1), 
p = -2 


10 


u r u t 


c — u, 
p = -2 






c = 0, 


1 1 
1 1 


-4(p+l)/(p-2) P ~~ 2 r „. \_. 

v p + 2 ; 


p^-2 






c = 0, 


12 


r~ 2 (3ru r — u){tu t — w) 


4fc = 6(1 -n)(p- 2), 
p = -4 



Table 10. Generalized-energy conservation laws for (n— l)(6 2 (p+l) 2 + c 2 ) 7^ 
0, p((6 + a) 2 + 6 2 c 2 (p + l) 2 ) ^0. 

in the next subsection 14.2.11 In contrast to the situation for the radial wave equation (l5a|) . 
some of the symmetries here are of contact type. Interestingly, all of the generalized-energy 
multipliers corresponding to the contact symmetries have a product form distinguished by 
the relations b = 0, or p = —2, or m = 4(p + l)/(2 — p). We will show later in subsection 
14.2.21 that these multipliers for the radial wave equation ( 15b)) can be derived alternatively 
from a variational type of Morawetz radial identity which produces conservation laws from 
symmetries in a similar way to the generalized-energy conservation laws that arise in the 
non- variational case of the radial wave equation fl5a|) discussed in subsection ( |4.1.2p . 
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4.2.1. Noether correspondence for u u — (c + b{p + l)u^)u Tr + {n — l)(cu r + bu^. +1 )/r and 
Uu = ku^u rr + (n — l)bu^ +1 /r. In the two variational cases, k = (p+l)b and c = 0, multipliers 
for conservation laws of the radial wave equation (15bj) are in one-to-one correspondence to 
variational symmetries of the respective Lagrangians (llObj) and (TT6T) . These cases can be 
merged if we consider the Lagrangian functional 

rti r<x> r 

C= / Lr m drdt, L = \(-u 2 t + cu 2 r ) + F(r,u r )du r (134) 

Jto Jo J 

given in terms of the expression (j!5p for F, with the coefficients b, c, k satisfying the condition 

c(k- (p+l)6) = 0. (135) 

Then Noether's theorem states that 

X = Pd/du (136) 

is a variational symmetry satisfying 

r m prXL = D t i4 + D r B (137) 
for some scalar functions A and -B depending on t, r, u, u t , and r-derivatives of u, u t , iff 

= - gu rr - f)r m P = D t T + D r X (138) 

is the characteristic form for a conservation law given by the conserved density and flux 

T = A- Q^—, X = B-Q^- (139) 

where 

Q = r m P (140) 

is the multiplier. 

For the generalized-energy multipliers in Table [TDJ we see that either Q is linear in ut and 

Q = u t a (t, r) + u r f3 1 (t, r) + uf3 (t, r) (141) 

where ao, @i, 0o have no dependence on u, in which case the corresponding variational 
symmetry is a restricted type of point symmetry, or Q involves a product of u t and u r , 

Q = u t (u r ax(t,r) + a (t,r,u)) + u r px{t,r,u) + f3 (t,r,u) (142) 

where a.\ has no dependence on u, while a and f3\ are at most linear in u, and f3 is at 
most quadratic in u, in which case the corresponding variational symmetry X is a restricted 
type of contact symmetry. This correspondence is summarized in Tables [TT] and [TJJ where 
C denotes the conserved quantity (|25p given by the conservation laws, and X denotes the 
variational symmetry generator given by the canonical form 

X = r- m (-a d/dt - fad/dr + (3 d/du) (143) 
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X 


C 


conditions 


1 & 3 


-d/dt 


ra(v* + cu 2 r ) + Kn ~ ^/iKOK" 1 dr 
Jo ' m 


c(k - (p + 1)6) = 


2 


— td/dt — rd/dr 
— ud/du 


POO 

/ [\t{u 2 + cu 2 + 2bh(u r )) + (ru r — u)u t )r~ 2 dr 
Jo 


k = {p+l)b,c^ 
n = — 1 


4 


— (p + A)rd/dr 
+ (2m — p)ud/du 


r°° 9h(n — 1 1 
/ l((p(m + 3)+A)t(u 2 + Zb[n l) h(u r )) 
Jo m 

+ 2y(p + 4)ru r + (2m — p)u)Ut)r dr 


c = 
-4 


5 


-2td/dt - w<9/<9« 


POO 

/ \{3{tu t - u)u t + ktu~ 2 )r m dr 

JO 


c = 0, 
p = -4 


6 


—t 2 d/dt — tud/du 


poo 

/ \{3{tu t -uf + kt 2 u; 2 )r m dr 
Jo 


c = 
p = -4 


7 


—r~ m d/dr 


POO 

/ u r u t dr 


c — u, 
p = -2 








c = 0, 


Q 
O 


—2rd/ dr — ud / du 


POO 

/ Ut(2ru r — ujr~ 2 dr 
Jo 


k — ouyn — l ) / z, 
p = -4 


9 


td/dt — r hird/dr 
+ ud/du 


/'OO 

/ ~((2u — tut)ut + 2(rkir)u r ut 
Jo 

+ 2kt\nu r ^r~ l dr 


c = 0, fc = 6(n-l), 
p = -2 



Table 11. Generalized-energy conserved quantities and variational point 



symmetries for (n — l)(b 2 (p + l) 2 + c 2 ) ^ 0, p(/c 2 + b 2 c 2 (p + l) 2 ) 7^ 0. 



in the point case (I14ip . or by 

X = r~ m (^ - (u r ai + a )d/dt - (u t ai + p x )d/dr + (3 d/du 

&r + u,uA w + &r> + u,{ w + + u '1m + -w) dldUt 

V du dr du dr dr du dr J J 

(144) 

in the contact case ()142p . 

Conserved quantity C\(= C3) is an energy arising from the time-translation symmetry 
Xi(= X 3 ). Conserved quantities C 2 , C 4 and C 5 are similarity energies arising from the 
scaling symmetries X 2 , X 4 and X 5 , while conserved quantity C 9 is a logarithmic counterpart 
of C4. Conserved quantity Cg and C-j are radial momenta arising respectively from the radial 
scaling symmetry X§ and the non-rigid dilation symmetry X7. Conserved quantity Cq is a 
conformal energy arising from the temporal inversion symmetry 

We will now state the vector space structure spanned by these conserved quantities as well 
as the algebra structure determined by the corresponding variational symmetries. 
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Theorem 9. For a multi-dimensional, hyperbolic, nonlinear radial wave equation (15b j) that 
has a variational structure (Q when a = pb ^ 0, c ^ 0, and an alternative variational 
structure (|T6j) when c = 0, variational point symmetries are admitted in the following cases, 
corresponding to all generalized-energy conservation laws whose multipliers are at most linear 
in derivatives. 

(i) For c = 0, the variational point symmetries comprise the time-translation X 3 ; the scaling 
symmetry X 4 ; the temporal scaling and inversion symmetries X 5 and X 6 when p = —4; the 
radial scaling symmetry X 8 when p = —A, k = 3b(n — l)/2; the non-rigid dilation symmetries 
X 7 whenp = —2 andXg whenp = —2, k = b{n—l). In the subcase p 7^ —2, p 7^ —A, the time- 
translation X3 and the scaling X4 generate a two-dimensional algebra with the commutator 
structure 

[X3,X 4 ] = -(p(^ + 3)+4)X 3 . (145) 

The corresponding two-dimensional vector space of conserved quantities is spanned by the 
energy C 3 and the similarity energy C 4 . In the subcase p = —4, k 7^ 3b{n — l)/2, the scaling 
X5 and temporal inversion X$ enlarge the commutator structure f!145j) by 

[X 3 , Xb] = -2X3, [X 4 , Xb] = 0, [X 3 , X 6 ] = X 5 , [X 4 , X 6 ] = 4(m + 2)X 6 (146) 

generating a four- dimensional algebra whose corresponding vector space of conserved quan- 
tities is spanned by C 3; C 4; and the similarity energy C5 in addition to the conformal energy 
Cq. In the subcase p = —4, k = 3b(n — l)/2, the radial scaling X 8 enlarges the commutator 
structure f)145p and f)146p by 

[X 3 , X 8 ] = [X 4 , Xg] = [X 5 , Xg] = [X 6 , X 8 ] = (147) 

generating a five- dimensional algebra whose corresponding vector space of conserved quanti- 
ties is spanned byC^, C 4 , C$, C%, andC%. In the subcase p = —2, k 7^ b(n — 1), the non-rigid 
dilation X 7 enlarges the commutator structure f)145|) by 

[X 3 , X 7 ] = 0, [X 4 , X 7 ] = 2(m + 1)X 7 (148) 

generating another three-dimensional algebra, whose corresponding vector space of conserved 
quantities is spanned by C 3 , C 4 and the radial momentum C 7 . In the subcase p = —2, 
k = b(n — 1), the non-rigid dilation symmetry Xg enlarges the commutator structure (11451) 
and (TT48D by 

[X 3 , Xfl] = X 3 , [X 4 , Xb] = -2X 7 , [X 7 , X 9 ] = -X 7 (149) 

generating a five- dimensional algebra whose corresponding vector space of conserved quanti- 
ties is spanned by C 3; C 4; C 7; and the logarithmic similarity energy Cg. 

(ii) For c 7^ 0, k = (p + 1)6 7^ ; the variational point symmetries comprise the time- 
translation Xi, and the scaling X2 when n = — 1. In the subcase n 7^ —1, the time-translation 
X 4 generates a one- dimensional algebra, while the corresponding vector space of conserved 
quantities is spanned by the energy C\. In the subcase n = —1, the scaling X2 and the 
time-translation Xi generate a two-dimensional algebra with commutator structure 

[X 1 ,X 2 ] = -X 1 (150) 

while the corresponding vector space of conserved quantities is spanned by C\ and the simi- 
larity energy C 2 . 
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X 


C 


conditions 




-r~ m (u r d/dt + u t d/dr 


POO 


c = 0, 
p= -2 


10 


+ u r u t d/du 


/ {\u 2 — k{l — \n.u r ))u r — b(n — l)u/r dr 
Jo 




+ mr~ x u r u t djdu^) 




11 


{{p - 2)u r - [p + 2)u)d / dt 
+ (p — 2)u t d/dr 
+ (p — 2)u r u t d / du 


— r(2k [ h(u T .)du r + u r u 2 )]r m dr 

p + 2 v J v ; tJ J 


c = 
4A; = 


-6(p-2)(n-l), 




+ (p + 2)u 2 d/du t 


p ^ - 


■2 




+ + 2)u r Ufd 1 1 du r 








(m — 3ru r )d/dt 








12 


+ 3r(u — tu t )d/dr 
+ (w 2 — 3rtu r u t )d / du 


POO 

/ j (2(2u — tu t ) (u — 3ru r )u t 
Jo 


c = 
4A; = 


-6(p-2)(n-l), 




+ (u — tu t )u t d/du t 


+ bt(n — l)(9ru r — u)u~ 2 )r~ 2 dr 


P = - 


■4 




+ (2tut — 3ru r — u)u r d / du T 









Table 12. Generalized-energy conserved quantities and variational contact 
symmetries for (n — l)b(p + 1) ^ 0, pk ^ 0. 



We will now state the vector space structure spanned by these generalized-energy con- 
served quantities as well as the algebra structure determined by the corresponding variational 
contact symmetries. 

Theorem 10. For a multi-dimensional, hyperbolic, nonlinear radial wave equation (15bj) that 
has a variational structure ffl6|) when c = 0, variational contact symmetries are admitted in 
the following cases, corresponding to all generalized-energy conservation laws whose multi- 
pliers are nonlinear in derivatives. 

(i) Forp = —2, the variational contact symmetries compriseXio generating a one- dimensional 
algebra, while the corresponding vector space of conserved quantities is spanned by the generalized- 
energy C 10 . 

(ii) Forp 7^ —2, Ak = — b{n — l)(p — 2), the variational contact symmetries comprise X n and 
Xi2- In the subcase p ^ —A, the symmetry Xn generates a one- dimensional algebra, while 
the corresponding vector space of conserved quantities is spanned by the generalized-energy 
C\\. In the subcase p = —A, the symmetries X12 and Xn generate a two-dimensional algebra 
with the commutator structure 

[X n ,X 12 ] = (151) 

while the corresponding vector space of conserved quantities is spanned by generalized- energies 
C12 and Cxi. 

4.2.2. Morawetz-type variational identities for Uu = (c + ku^)u rr + (n — l)(cu r + bu^. +1 )/r. 
We will now derive two different Morawetz identities of variational type for the radial wave 
equation ( l5b]) in the case c = with k 7^ 0, b(p + 1) 7^ 0. The steps are similar to the 
derivation of the identity (11191) for the non- variational case of the radial wave equation (JJa]), 
and so we omit most details. 
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To derive the first identity, we use the Lagrangian ( TT6|) to write 

r m (u tt - ku p r u rr - (n - l)bu p r +1 /r) = E u (r m L ) + E u {r m L x ) (152) 

where E u is the 1st order Euler-Lagrange operator (11001) (i.e. local variational derivative) 
with respect to u, and where 

k 

L = -\u\, Li = -^-j—^ h ( u r) (153) 

are the temporal and radial parts of the Lagrangian L = L + L\. We now multiply equation 
(11521) by the expression 

W = ru r - Xu (154) 

which represents the characteristic function given by the generator of a scaling transformation 
on (r, u), namely St = 0, Sr = r, Su = Xu, for some constant A. By using the operator identity 
( 11021) combined with the radial scaling properties of Lq and L±, we obtain 

(ru r - Xu)E u (r m L ) = -E u (r m L ) + r m (2rD r + (m - A + 1))L (155) 

and 

(ru r - Au)S u (r m Li) = -S u (r m Li) + r m (p - m + 1 - A(p + l))Li (156) 

with 

Lo = |(ra r — \u)v,t (157) 

and 

~ k [ ~ k ~ 

L\ = (2r / h(u r )du r +(Xu — ru r )h(u r )) + 9ru r = (ru r h(u r ) + Xuh(u r )) + (9 + 2)ru r 

p + 1 J p+3 

where 9 = (A — l)/(m + 1) when p = —2 and = otherwise. Finally, we require that the 
terms in equation (11561) containing L\ vanish. This yields 

A = 1 - ( n - l)b/k. (159) 

Thus, we have expressed the Euler-Lagrange terms (I155P and (1 1 5 6 [) as the sum of a total 
variational derivative E u (r m L) for a new Lagrangian 

L = L + L h (160) 

plus a dilational derivative (rD r + X)(r m L /2) with 

A = (1 — m - A)/2 = b(n - 1)(1 - p/2)/k. (161) 

As a result, equation fl 1 5 2 [) multiplied by expression (I154p yields the variational dilation 
identity 

r m {ru r - Am) (u u - ku p r u rr - (n - l)bu p+2 /r) = -E u {r m L) - (rD r + X)(r m u 2 t ). (162) 

This identity fl 1 6 2 [) is able to produce conservation laws from multipliers that are given by 
certain symmetries X = Pd/du in characteristic form, in a manner similar to Noether's 
theorem (l39p -( l4"Uj) . In particular we consider the time translation 

X = utd/du. (163) 

When the identity ( I162p is multiplied by 

P = u t (164) 
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we see first that the Euler-Lagrange term yields a total time derivative 



fjj _ AT 

- u t E u (r m L) = D t (r m u t —) - r m W XL = D t {r m u t — - r m L) (165) 

since, by inspection, the Lagrangian functional £ = J* t 1 J °° Lr m drdt is invariant under time- 
translation. Second, we see that the dilational derivative term produces 

-u t {rD r + X){r m u 2 ) = -2r m+1 u 2 t u tr - (A + m)r rn u 3 t = D r ( - §r m+ V) + (1(2 - m) - ~X)r m u 3 t 

(166) 

which will reduce to a total radial derivative iff 

A = (2-m)/3. (167) 
Equating expressions ( I167p and (11611) . we obtain the condition 

Ak = -b(n - l)(p - 2) (168) 

and hence 

A = -(m + l)/3, m=-4(p + l)(p-2). (169) 

Thus, for a, b, k, A satisfying conditions fl 168[) and (I169p . the product of the variational 
dilation identity (I162p and the symmetry characteristic function (I164p yields a conservation 
law, where the conserved density and flux are given by T\\ and X\\ in Table [9] while the 
multiplier is given by Qu in Table [TUJ We note that the form of this multiplier consists of 
r m times the product of the characteristic functions W = ru r — Xu and P = u t representing 
the scaling transformation and the time-translation symmetry respectively. The similar 
conservation law given by the conserved density T 12 and flux X 12 in Table [9J with the 
multiplier given by Q12 in Table [TOl can be obtained from the variational identity (I162p by 
using the temporal scaling symmetry with the characteristic function P = tu t — u (in place 
of the previous time-translation symmetry). 

To derive the second identity, we omit the radial factor r m from equation (I152p . which 
then yields 

u tt - ku p r u rr - (n - l)K +1 /r = E U (L ) + E U {L X ) + /i< +1 /r (170) 

where 

/i = -(n-l)6 (171) 

is the coefficient of the non-variational term. Now, instead of multiplying this equation by 
the scaling characteristic function (I154p . we use 

W = u r (172) 

which represents the characteristic function given by the generator of a translation on r. We 
thereby obtain 

u r E u {L ) = -E U {L ) + D r L (173) 

and 

u r E u {L x ) = -£ u (Li) (174) 

with 



k .. . { k 
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\u r u 2 , Li = (2 / h(u r )du r — u r h(u r )) = u r h(u r ) + 9ru r (175) 

p+1 J p+3 



where 9 = Ik when p = —2 and 8 = otherwise. Substituting expressions (I173p and (I174p 
into equation (I170p multiplied by the characteristic function (j!72p . we find 

u r (u tt - ku p r u rr - (n - l)bu P . +2 /r) = -E U {L) - D r (\u 2 t ) + /iu p r +2 /r. (176) 

Finally, we observe that the term u p+2 /r in equation (I176p will be a total radial derivative 
iff 

p = -2 (177) 

in which case 

HU*+ 2 /r = D T (jilnr). (178) 
As a result, equation (I176P reduces to the variational identity 

u r (u tt - ku p T u rr - (n - l)bu p+2 /r) = -E U {L) - D r (\u 2 t + (n - 1)6 In r) (179) 

where 

l = U + l x (180) 

is a new Lagrangian. This identity (11791) is able to produce conservation laws in a similar 
way to the previous identity (I162p . In particular, we consider the time translation (I163p . 
Multiplying the identity (I179P by the characteristic function (11641) . we see first that the 
Euler-Lagrange term yields a total time derivative 

dL ~ BL 

- u t E u (L) = D t (u t —) - prXL = D t (u t — - L) (181) 

since, by inspection, the Lagrangian functional £ = J °° Ldrdt is invariant under time- 
translation. Second, we see that the radial derivative term produces the total derivatives 

- u t D r {\u 2 t + (n - l)Hnr) = D r (- ^uf) + D t (-(n - l)bu/r). (182) 

Thus, for p = —2, the product of the variational identity (I179P and the symmetry char- 
acteristic function (I164p yields a conservation law, where the conserved density and flux are 
given by Tio and Xio in Table |9] while the multiplier is given by Qio in Table [lOj We note 
that the form of this multiplier consists of just the product of the characteristic functions 
W = u r and P = Ut- Interestingly, since the Lagrangian (I180p is explicitly invariant under 
the shift symmetry X = d/du with the characteristic function P = 1, we can use the vari- 
ational identity (11791) to obtain the conservation law given by T 7 and X 7 given in Table El 
whose multiplier is shown in Table [TUJ 

4.3. Computational remarks. 

In the computation of conservation laws for a radial wave equation ([3]), the multiplier 
determining system (l30l)-(l3~ll is formulated in the jet space whose coordinates are defined by 
t, r, u, u t , u r , u tr and u rr , where u u is replaced by gu rr + f. For both wave equations ( )5a| and 
( l5b|) . this determining system splits with respect to u t , u tr and u rr into an overdetermined 
system of 13 equations for unknowns Q = a(t,r,u,u r )u t + (3(t,r,u,u r ), a, b, c, p and n. 
The resulting overdetermined systems are linear in the unknown multiplier components a 
and /3 but are nonlinear jointly in all the unknowns. In particular, the unknown p appears 
non-polynomially in each overdetermined system through respective exponents of u and u T . 

We again carry out all computations by means of computer algebra. The program Con- 
LAW4 [T5] is used to formulate and split the conservation law determining system to get the 
overdetermined system which in turn is solved for the 7 unknowns a, (3, a, b, c, p and n by 
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the package Crack [19]. For both wave equations (l5a| and (|5b|) . Crack is able to complete 
the entire computation of multipliers automatically. Compared to the computation of point 
symmetries (cf. section [373]) . the computational difficulties are similar but the total number 
of steps needed is approximately 2 to 4 times greater. 

Once all solutions of the overdetermined system for multipliers Q are obtained from 
Crack, we compute the conserved densities T and fluxes X interactively in Maple 14 
by evaluating the homotopy integrals (I34p -( 13T)]) whenever possible or by otherwise solving 
the linear PDE system (132]) - (133]) . All conditions on the parameters a, b, c, p and n for a 
multiplier to exist are imposed at the start of the Maple computation. 



In this paper we have fully classified all point symmetries ( I19p . kinematic conservation 
laws (l56]) - (l57j) . and generalized-energy conservation laws (|46]) - (|47]) for the two classes of 
nonlinear radial wave equations (15a]) and (15bj) which are parameterized in terms of the 
constant coefficients a, b, c, and the nonlinearity power p, as well as the dimension n. The 
only restrictions placed on these five parameters are that we want each wave equation to be 
multi-dimensional, hyperbolic, and semilinear or quasilinear. 

Since Noether's theorem is applicable for the cases in which these wave equations have a 
variational formulation, our classification of conservation laws directly yields a corresponding 
classification of all variational point symmetries and also variational contact symmetries 
of a certain restricted form. We find that such variational cases account for all of the 
generalized-energy conservation laws admitted by radial wave equation (15bl) but only some 
of the generalized-energy conservation laws admitted by radial wave equation (l5aj) . 

In the non- variational cases for radial wave equation (|5a]) . we show that the additional 
generalized-energy conservation laws arise through a new type of radial dilation identity 
that produces conservation laws from symmetries in a different way than Noether's theo- 
rem. A similar radial dilation identity for radial wave equation (15b j) is able to reproduce 
the generalized-energy conservation laws that correspond to all of the variational contact 
symmetries. 

Some of the symmetries and conservation laws found in our classification involve non- 
integer values of n. We can interpret such cases by rewriting the wave equations (l5aj) and 
(I5bj) in a 2-dimensional form in terms of a parameter v = 2 — n as follows: 



5. Concluding Remarks 





and 





where 





and 




(185) 
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This form of the wave equations (!5a| and (15bl) is applicable for any value of n. The 
meaning of v can be understood through the kinematic quantity 



C{t) = I urdr (186) 

which satisfies 



o 



d 2 C 
dt 2 



uN + S (187) 



if lim u r = 0, where 



poo 

N = limYcw + bH(u)), S= (a - pb)u p ~ l ulrdr (188a) 

Jo 

for equation (I183al) . and 

POO 

N = \imcu, S= b(l - v)u p r +1 - kH(u r )dr (188b) 

Jo 



for equation (1183b|) . We see from equation (I187P that if u(t, r) is viewed as the amplitude of 
a vibrating surface then C(t) describes the net transverse displacement of the surface and 
hence vN has the interpretation of a time- dependent forcing term applied to the surface at 
r = 0. 

The symmetries and conservation laws we have obtained in Tables [T] to [12] in sections [3] 
and H] can be expected to have several important applications in the analysis of solutions 
for the n-dimensional radial wave equations (jS]) (in the case of positive integer values of 
n) or the radial surface wave equations (11831) (in the case of all other values of n). For 
example, in variational cases, both of the n-dimensional radial wave equations ([5]) possess a 
conserved energy as well as a conserved similarity energy for c = and a conserved conformal 
energy for c = 0, p = —4. The multipliers for these energies can be used to derive useful 
energy identities in all non- variational cases for these radial wave equations. Other typical 
applications for our results have been mentioned in section [1] and will be explained in detail 
elsewhere. 

There are several directions in which the work in this paper can be extended. First, we 
plan to generalize our classification of symmetries and conservation laws to include all contact 
symmetries (T2"TT) and all 1st order conservation laws (12TJ1) for the radial wave equations §5§ in 
all cases n^l. Second, in a sequel to the present paper, we plan to determine which of the 
symmetries and conservation laws for the n-dimensional case of the radial wave equations ([5]) 
can be lifted to yield corresponding symmetries and conservation laws for the translationally- 
invariant n-dimensional wave equations (jl]). 
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